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Abstract 

The two-by-two Sp(2) matrix has three parameters with unit de- 
terminant. Yet, there are no established procedures for diagonalizing 
this matrix. It is shown that this matrix can be written as a similar- 
ity transformation of the two-by-two Wigner matrix, derivable from 
Wigner's little group which dictates the internal space-time symme- 
tries of relativistic particles. The Wigner matrix can be diagonalized 
for massive and space-like particles, while it takes a triangular form 
with unit diagonal elements for light-like particles. The most imme- 
diate physical application can be made to repeated one-dimensional 
transfer matrices appearing in many different branches of physics. An- 
other application of current interest could be the dis-entanglement of 
entangled systems. 
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1 Introduction 



The two-by-two matrices of the group Sp(2) or its unitary equivalent appear 
in almost all branches of physics. Indeed, a systematic approach to the study 
of those matrices is a respectable branch of physics [U [2]. 

The Sp(2) matrix is the basic language for linear canonical transforma- 
tions for phase-space approach to both classical and quantum mechanics [3]. 
The same is true for classical and quantum optics [3], as well as the current 
problem of entangled harmonic oscillators [H [5] . The one-dimensional scat- 
tering matrix also takes the form of this Sp(2) matrix. The two-by-two beam 
transfer matrix is unitarily equivalent to that of the Sp{2) group 0, El IE]- 
The list could be endless. In short, we cannot do modern physics without 
two-by-two matrices, especially without those of the Sp{2) group. 

Do we then know how to diagonalize this simple mathematical expression? 
We can if they can be diagonalized by a rotation. If not, there are no standard 
procedures for approaching this problem. In this paper, we address this 
fundamental issue. Our main interest is whether this matrix can be written as 
a similarity transformation of a diagonal matrix. If that is the case, we solve 
the practical problem of calculating repeated applications of the same matrix 
as we see in finite one- dimensional crystals, including periodic potentials [61 
IE], multilayer optics [HIHD], lens optics [HUE], an d laser cavities [T3], and 
many other problems in physics. 

We first show that the most general form of the Sp(2) matrix can be 
written as a similarity transformation of Wigner's little-group matrix [14] . 
This Wigner matrix takes three different forms depending on the parameters 
of the original Sp(2) matrix. However, one of these three matrices is not 
diagonalizable, but offers the same conveniences as diagonal matrices do. 
Indeed, the diagonalization of the Sp(2) matrix requires the construction of 
a two-by-two representation of Wigner's little groups. 

For this purpose, we start with the Bargmann decomposition of the Sp(2) 
matrix in which it is written as a product of three one-parameter matrices. 
We also write the matrix as a similarity transformation of Wigner's little 
group matrix. In this way, the Wigner parameters can be written as the 
Bargmann parameters, allowing us to write the Sp(2) matrix as a similarity 
transformation of the Wigner matrix, namely Wigner's little-group matrix. 

In section [2], we start with the Bargmann representation and transform 
it into a form convenient for writing it as a similarity transformation of 
the Wigner's little-group matrix discussed in section [3J In section 0J the 
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Bargmann parameters are written in terms of the Wigner parameters. The 
Sp{2) matrix takes many different forms. They are shown to be unitarily 
equivalent in section [5j In section [6], we discuss some immediate applications 
of the diagonalization by Wignerization of the Sp(2) matrix. 



2 Bargmann decomposition 

The Sp(2) matrix appears in the physics literature in several different forms. 
This two-by-two unimodular matrix is often written as 

M=d n), (1) 



C D, 

with (AD — BC) = 1. All four elements are real numbers. This matrix is 
commonly called the ABCD matrix and has three independent parameters. 
Furthermore, this unimodular matrix M can also be written as 

M = R 1 (6 1 )B(2X)R 2 (6 2 ), (2) 

with 

/ cosh A sinh A 
V sinh A cosh A 



B(2X) 



/cos(^/2) -sin(^/2)\ 

UiM^/2) cos^/2) J' [6) 

This is known as the Bargmann decomposition [15J. We use the notation 
B(2X) instead of B(X) purely for convenience. The main point is that the 
original three-parameter matrix is decomposed into three one-parameter ma- 
trices. 

We can rewrite the Bargmann decomposition of equation (T5]) as 

M = R 1 BR 2 = (LR) B fa- 1 ) = L(RBR)L-\ (4) 

where 

#i = LR, R 2 = RL~\ (5) 

and 
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with 

6>i + #2 . Oi — 9 2 

9 = , o = . 

2 ' 2 

Then, the matrix M can be written as a similarity transformation of (RBR) 
with respect to L, or a rotation of the system by L. Thus, for all practical 
purposes, it is sufficient to study the core matrix RBR with two parameters. 
This is still a form of the Bargmann decomposition. 

Now the question is how to diagonalize the core matrix (RBR) and write 
it as a similarity transformation. This core matrix takes the form 

R{9)B{2\)R{9) = ( ( [ CO * hX) r -(coshA)sin^ + sinhA\ 
V (cosh A) sin 9 + smh A (cosh A) cos 6 J 

We assume here that both sm9 and sinhA are positive. 

The product of the two off-diagonal components becomes 

cosh 2 A[tanh A + sin 6} [tanh A - sin 6} = (cosh A) 2 cos 2 9-1, (8) 

Thus, if one of the off-diagonal components vanishes, the diagonal elements 
become one, and the RBR matrix becomes 

ff + ( 8tahA '=(2smhA !)■ 01 (J ~ 2 T hA ). 0) 

respectively. This form is known as the Iwasawa decomposition, and is a spe- 
cial case of the Bargmann decomposition. This matrix is not diagonalizable. 

Yet, we are interested in diagonalizing the RBR matrix by following the 
usual procedure of calculating the eigenvalues of the matrix, which leads to 
the eigenvalues 



E± = (cosh A) cos #± ^/(coshA) 2 cos 2 # - 1. (10) 

with E + E_ = 1. The quantity inside the square root sign is the same as 
the product of the two off-diagonal elements given in equation (jSJ) . Then the 
diagonal matrix takes the form 

If (cosh A) 2 cos 2 9 is smaller than one, and the quantity inside the square- 
root sign is negative, the eigenvalues are 

E± = exp (±i0/2), (12) 
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with 

cos(0/2) = (coshA)cosfl. (13) 
If (cosh A) 2 cos 2 9 is greater than one, the eigenvalues become 

£ ± = exp(± X /2), (14) 

with 

cosh(x/2) = (cosh A) cos^. (15) 

If (cosh A) cos 9 = 1, the quantity inside the square-root sign vanishes, 
the eigenvalues collapse to one, and the RBR matrix becomes one of the 
triangular matrices of equation (jUJ). 

In view of the triangular matrices of equation (Q, the best we can do is 
to write the RBR matrix as a similarity transformation of 

( elH2 ^ ( eX/2 \ r™\ 
{ o e -^ 2 J' I e-W' { ' 

or one of the triangular matrices of equation We shall call the set of 
these three matrices the Wigner matrix, use the notation W± for the set of 
three matrices with H±, and W collectively for both. 

Our plan is two write RBR as a similarity transformation of the Wigner 
matrix: 

RBR = S^WS. (17) 

As we shall see in Sec. [31 the W matrices are also derivable from Wigner's 
little group which dictate the internal space-time symmetries of elementary 
particles [H] 

Then the question is which among the two Wigner matrices is to be 
chosen. In order to address this question, let us note the Lorentz group has 
two branches. The Sp(2) group is generated by 

forming a closed set of commutation relations: 

[J 2 , K 3 ] = iK u [J 2 , K x \ = -iK 2 , [J 2 , K 3 ] = -U 2 . (19) 

The rotation generator J 2 is antisymmetric, and K 3 and K\ are symmetric. 
Furthermore, this set of commutation relations is invariant under the sign 
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change of the K~i matrices while it is not when J 2 changes its sign. The rota- 
tion matrices generated by J 2 matrix are anti-symmetric, while the squeeze 
matrices generated by K~i are symmetric. 

Let us go back to the RBR matrix of equation ([7j). Since the B(2\) 
matrix is symmetric and is generated by K\. We should therefore consider 
the case where B{2\) is replaced by B(— 2A). The resulting RBR matrix is 

R(9)B( 2X)R(9) — ( (cosh A) cos # — (cosh A) sin6> — sinh A 

V (cosh A) sin^ — sinh A (cosh A) cos 9 

(20) 

For convenience, we shall use the notations (RBR) + and (RBR)- for equa- 
tions (j7]) and ([20]) respectively: 

(RBR)+ = R(9)B(2\)R(0), 

(RBR)_ = R(9)B(—2X)R(9). (21) 

We shall still use the notation RBR collectively for (RBR) + and (RBR)-. 
We cannot obtain (RBR)_ from (RBR) + by simply changing the sign of the 
A parameter. It is a parity operation. 

If both sin# and sinh A are positive, the upper-right element of (RBR) + 
can vanish, while the lower- left element can become zero for (RBR)-. They 
are therefore consistent with i/ + (sinhA) and if_(sinhA) respectively. 

Let us see why the W matrix should be called the Wigner matrix in 
section [31 



3 Wigner's little groups 

Since the group Sp(2) is locally isomorphic to the Lorentz group applicable 
to two space dimensions and one time variable, we can import useful ma- 
trix identities from the kinematics of the Lorentz group, particularly from 
Wigner's little group [H]. 

Wigner's little group is the maximal subgroup of the Lorentz group whose 
transformations leave the given four-momentum of a particle invariant |14j . 
Let us first consider a particle moving along the negative z direction with 
the four-momentum 

P = (E,0,0,-P), (22) 
with E = y/p' 2 + m 2 . We use here the four- vector convention (t,x,y,z), and 
the c = H = 1 unit system. 
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It is possible to bring this four-momentum to 



(m, 0,0,0), 



(23) 



by applying the transformation matrix 



( cosh 77 



V sinh 77 




1 







1 




sinh 7] \ 



cosh r\ J 



(24) 



width 



E = m(coshr/), p = m(sinhr/). 



(25) 



The four-momentum of equation ([23]) is invariant under the rotation 



R(<f>) 



(I 



Vo 





cos (f) 


— sin . 



\ 

sin (f. 

1 
cos 4> ) 



(26) 



After this rotation, we can bring back the four momentum of equation (|23|) 
to the starting momentum of equation (1221) by applying the inverse of the 
boost matrix S(rj). 

Thus the transformation S(—r])R(<j))S(r]) will leave the four-momentum 
of equation (122]) invariant. 

We can also rotate the entire system around the z axis without chang- 
ing the four- momenta of equation (12"!?]) and equation f[2"3"]) . and the matrix 
S(r]). Thus, we can obtain the three-dimensional rotation group by tak- 
ing into account this degree of freedom along with the rotation matrix of 
equation (|26j) [15] . The group represented by this Lorentz-boosted rotation 
matrix is known as Wigner's 0(3)-like little group for massive particles. 

If the four-momentum is space-like, it can be brought to (0, 0, 0, p). 
This four-momentum is invariant under the boost 



x{x) 



( cosh x 
sinh x 


V 



sinh x 
coshx 










1 








1/ 



(27) 



along the x direction. The little group is thus represented by S~ 1 X(x)S, 
along with rotations around the z axis. Wigner's little group for space-like 
four-momentum is a Lorentz-boosted 0(2, 1) group. 
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If the particle is massless, there are no Lorentz frames in which the particle 
is at rest, but its four-momentum can be brought to (k, 0, 0, k). This 
four-vector is invariant under rotations around the z axis. In addition, it is 
invariant under 



iV( 7 ) 



/l + 7 2 /2 -7 

- 7 1 



7 2 /2 



7 





1 




-7 2 /2 \ 

7 


l- 7 2 /2/ 



(28) 



In his 1939 paper, Wigner showed that this N operator together with the 
above-mentioned rotation around the z axis form a three-parameter group 
which is locally isomorphic to the two-dimensional Euclidean group. Its 
rotational degree of freedom corresponds to the helicity of the particle while 
the two-translational degrees collapse into one gauge degree of freedom [T7] . 

This summarizes what Wigner did in his 1939 paper and later papers 
on massless particles [HI [17]. The transformation matrix which leaves the 
four-momentum invariant is 



S(-r,)WS(r,), 



(29) 



where W takes the form of R(4>), X(x), or iV( 7 ). 

It is well-known that the six-parameter Lorentz group can also be repre- 
sented by two-by-two unimodular matrices. This group of two-by-two matri- 
ces is called SL(2, c), and Sp(2) is one of its subgroups. 

Wigner's approach is not the only method to obtain the momentum- 
preserving transformations. If we rotate the four-vector of equation (|22|) 
around the y axis using the rotation matrix 
















cos# 





sin 6? 








1 





\o 


— sin# 





cos 9 ) 



R{9) 



the resulting four-momentum is 

(E, —p sin 9, 0, — p cos ( 



(30) 



(31) 



The above four-by-four matrix corresponds to the two-by-two matrix R(9) 
of equation ([6]). 



8 



It is possible to boost this four- vector along the positive x direction using 
the boost matrix corresponding to B(2X) of equation 0. Then the four- 
momentum becomes 

(E,p sin 9,0, -pcos6). (32) 
The four-by-four boost matrix is 



B(2X) 



^cosh(2A) sinh(2A) 0\ 

sinh(2A) cosh(2A) 

10 

V 1/ 



(33) 



We can then rotate this four-momentum to the original form of equa- 
tion (12"2"|) . using again the rotation matrix R{9) of equation (13"U1) . This method 
was discussed in detail in Ref. [18], but only for a massive particle whose four- 
momentum takes the form of equation (T22j) . 

It should be noted that, in order that the four- vector of equation (1221) re- 
turn to the original form after the rotation-boost-rotation process, the boost 
parameter A cannot be arbitrary. It should be determined from the rotation 
angle 9 and the four-momentum of equation (122j) . 

On the other hand, if the parameters 9 and A are chosen first, they deter- 
mine the form of the four- vector. It could be time-like (non-zero massive), 
space-like (for imaginary mass), or light-like (zero mass), covering all three 
possible little groups [H] . 

We choose the notation W for the set of the four-by-four matrices R(<f>) of 
equation (126]) . X(x) of equation (127|) . and N(j) of equation (128]) . and call it 
the 'Wigner matrix." Then, in terms of W, the matrix RBR can be written 
as a similarity transformation 

R{9)B(2\)R{9) = S(-T])WS{r]). (34) 

We are then interested in whether the same relation can be derived for the 
two-by-two matrices of the Sp{2) group. 

Since the group Sp(2) is locally isomorphic to 0(2, 1), we expect to be 
able to write this similarity transformation in terms of the corresponding two- 
by-two matrices. Indeed, there are two-by-two Sp{2) matrices corresponding 
to the four-by-four matrices for R((p),X(x), an d N(j) given equation 
equation (T271) . and equation (T28]) respectively. They are 

/cos(0/2) - sin (0/ 2 )\ 
UW) Wn(0/2) cos(0/2) )> 
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x±{x) 



( cosh(x/2) ±sinh(x/2)\ 
V±sinh(x/2) cosh(x/2) ) 



NAl)=( 1 1 J), or N_(l)=(l x 7 ). (35) 
The two-by-two counterpart of the boost matrix of equation (|24"1) is 

^) = (T e-°/0" (36) 

Thus the two-by-two representation of Wigner's little group should also be 
written as S(^f7])WS(±rj). However, we have to confront the question of 
which N(^) matrix from equation (1351) is to be chosen. If we go back to 
the (RBR)± matrices of equations ([7]) and (|20|) . the upper-right element can 
vanish for (RBR) + , while (RBR)^ can have a vanishing lower- left element. 

We shall also see that the transition from X + (x) to X_(— x) cannot be 
achieved through an analytic continuation from positive x to negative x m 
equation (1351) . We shall assume that A and x have the same sign. 

Indeed, the two-by-two Wigner matrix has two branches. Thus N + and 
X + should be chosen for (RBR) + . We shall call this set of Wigner matrices 
W+, and write 

(RBR) + = S(-r))W + S(r)). (37) 

What then happens to iV_ and X_? 

In order to answer this question, let us start with the four-momentum 

P=(E,0,0,p), (38) 

instead of (E, 0, 0, —p) of equation (1221) . We can boost this four-momentum 
to E = (m, 0,0,0), using the inverse of S(r]) given in equation (1211) . This 
four-vector is invariant under rotations. We can then boost back to the 
starting form of equation (|38|) . 

As for the Bargmann decomposition, we can rotate the four-vector of 
equation (I3"8"l) to 

P = (E,p sm9, 0,p cos 9), (39) 

by applying the rotation matrix R(8). We then apply the boost matrix 
B(— 2A) to this four vector to get 

P = (E,-psin9,0,pcos0). (40) 
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The rotation matrix R(9) will bring this vector back to the original form of 
equation (j38p . If we complete this kinematics, the resulting matrix is (RBR)_ 
of equation fl20|) . This version of the Lorentz kinematics was thoroughly 
discussed in reference [T8] . 

The resulting RBR matrix is that of equation (1201) . which we write as 
(RBR)_. The lower-left element of this matrix can vanish. The correspond- 
ing Wigner matrix should include iV_ instead of N + , and 

(RBR)_ = S(rj)W-S(-Ti). (41) 

If we combine equation f[3"T|) and equation (l4"Tj) . we can write 

(RBR) ± = S(TV)W ± S(± V ). (42) 

4 The Sp(2) matrix as a similarity transfor- 
mation of the Wigner matrix 

Our main concern is to see whether the M matrix or its core RBR of equa- 
tion (j7]) or equation (1201) can be written as a similarity transformation of 
the diagonal D matrix of equation ffTTj) . We have already shown that this is 
not always possible when the matrix is triangular. On the other hand the 
triangular matrix has one parameter, and is as simple as the diagonal matrix. 

We propose to write this expression as a similarity transformation of the 
Wigner matrix W given in equation (1351) . Let us combine RBR matrices of 
equation (J7]) and equation (120]) into one expression: 



R{6)B{±2\)R{6) 



{ (cosh A) cos 9 —(cosh A) sin 9 ± sinh A ' 

V (cosh A) sin# ± sinh A (cosh A) cos 9 

(43) 

If the diagonal elements are smaller than one, then 

(cosh A) sin 9 > sinh A, (44) 

and the off-diagonal elements of equation (1431) have opposite signs, the RBR 
matrix can be written as 

R{9)B{±2\)R{9) = S{tv)R{<P)S{±v)- (45) 
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The right-hand side takes the form 

cos(0/2) -e^sin(0/2)\ 
e ±r > sin(0/2) cos(0/2) J ' 1 ; 



By comparing this expression with equation (]43|) . we can write and r\ in 
terms of A and 9, and the result is 

cos(0/2) = cosh A cos 9, 

2r] cosh A sin 9 + sinh A 
cosh A sin 9 — sinh A 



cosh(x/2) ±ef sinh(x/2) A 
±e ±J? sinh(x/2) cosh(x/2) )' 



If the diagonal elements of equation (1431) are greater than one, then 

(cosh A) sin 9 < sinh A, (48) 

and the off-diagonal elements have the same sign. We should then use X(x) 
as the Wigner matrix, and write RBR as 

R(9)B(±2X)R(9) = S(±r ] )X ± ( X )S(Tv)- (49) 
The right-hand side of this expression can be written as 

(50) 

If we compare this expression with the (RBR)± matrix of equation ( T4HD , 

cosh(x/2) = cosh A cos 9, 

9 „ cosh A sin 9 + sinh A 

e v = . (51) 

sinh A — cosh A sin 9 

The upper-right or the lower-left component of equation f|43l) can go 
through zero starting either from a small negative number or from a small 
positive number. 

For the first case, if (cosh A sin 9 — sinh A) goes through zero from a small 
negative to a small positive number, for the (RBR) + branch one should start 
from S(—r])R((j))S(ri). Then 

H + {2 sinh A) = S(-r])R((j))S(r]), (52) 
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hence in this negative region we have 



2sinhA = e''sra(0/2). (53) 
For the (RBR)_ branch one should start from S(r])X_(x)S(— rj), then 

#_(-2sinhA) = S(ti)X-(x)S(-ti). (54) 
So in this negative region we have 

2sinhA = e' 7 sinh(x/2). (55) 

When rj goes to infinity, the upper-right component of S(—r))R((f>)S(r]) 
section of equation ( 1461) . namely — exp~ v sin(0/2) goes to zero. Simi- 
larly, ± exp~ v sinh(x/2) vanishes, which is the lower-left component of 
S(t))X-(x)S(— rj) section of equation equation (j50|) . In those cases <p and 
X should be very small so that the diagonal elements collapse to one. This 
limiting process is known as group contraction, and is applicable to various 
problems in physics [T2j [T7]. 

On the other hand when (cosh A sin 9 — sinh A) goes through zero from a 
small positive number to a small negative number, on the (RBR) + branch 
one should start from S(—1])X+(x)S(t}), thus 

# + (2sinhA) = S(-r))X + (x)S(r)). (56) 

In this positive region we have the relation between the parameters as in 
equation (I55p . On the (RBR)_ branch we should start from S(r))R((j))S(—r)) 
thus 

#_(-2 sinh A) = S(rj)R(<f>)S(-7j). (57) 

So in this positive region we have the relation between the parameters as in 
equation (1531) . 

When (cosh A sin 9 — sinh A) = the diagonal elements are equal to one, 
one of the off-diagonal elements vanishes. The similarity transformation 
should be written as 



# ± (sinhA) = S(Tv)N±(j)S(±r]). (58) 

Thus 

2 sinh A = <ye v . (59) 
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Since there is only one variable for the Bargmann decomposition collapses 
into the Iwasawa decomposition, the variables 7 and 77 become combined into 
one variable. 

Let us now return to the original question of diagonalizing the ABCD 
matrix or its Bargmann decomposition given in equation (j2J). We now have 
RBR matrix as a similarity transformation of the W matrix. The transfor- 
mation matrix S is given above. Thus, the M matrix can be written as a 
similarity transformation given in equation (T4J). The transformation matrix 
is G = LS~ X , where the L matrix is defined in equation d^J). The G matrix 
now takes the form 



tqt - ( eTv/2 cos(5/2) -e ±r >/ 2 sin(5/2) \ 
~ ~ 1 e^/2 sin ( 5 / 2 ) e ±v/2 CO s(5/2) J ' 



(60) 



We thus conclude that the Sp(2) matrix cannot always be diagonalized, 
but it can be written as a similarity transformation of the Wigner matrix. 
The transformation matrix is given above. 



5 Further Mathematical Details 

The Sp(2) matrices are used in many different branches of physics. They 
take different forms in the literature, but they are unitarily equivalent. The 
most convenient way to organize those different expressions is to recognize 
that the group Sp(2) is a subgroup of the six-parameter SL(2, c) group, 
which is locally isomorphic to the Lorentz group applicable to three space 
and one time dimensions. The three-dimensional rotation group is one of its 
subgroups. In addition, it has three 0(2, 1) Lorentz groups applicable to two 
space and one time dimensions. 

The form given in equation (T5]) corresponds to the 0(2, 1) group applicable 
to z, x and t, where rotations around the y axis are allowed. On the other 
hand, the expression 

fe ie ^ 2 \ /cosh A sinh A \ / e id2 ' 2 \ ( 
{ e^/VUnliA cosh A A e"^/ 2 J {b } 

is quite common in the literature, especially in layer optics. This is also the 
form used by Bargmann in his original paper [15]. This expression represents 
the 0(2, 1) subgroup applicable to (x, y, t) while rotations around the z axis 
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are allowed. Indeed, the decomposition of equation (j2J) can be obtained from 
equation (1611) by a conjugate transformation with the matrix 

i n i 



V2U v (62) 

It is also possible to decompose the Sp(2) matrix in terms of the three 
matrices of the form 



cosfl -sin0\ fl - 7 \ fl 
sin^ cos# ) ' VO 1 ) ' V7 1 



(63) 



These expressions are useful in para-axial lens optics. They are discussed in 
detail in Ref. pj]. 

Let us go back to the ABCD matrix of equation ([1]). In terms of the 
Bargmann parameters, they can be written as 

A = — sinh A sin 5 + cosh A cos 9 
B = sinh A cos 5 — cosh A sin 9 
C = sinh A cos 5 + cosh A sin 9 

D = sinh A sin 5 + cosh A cos 9. (64) 
Conversely, the Bargmann parameters are 

n C-B . D-A 

tan6 = A + D > tan 6= bTc } ( ' 



and 



< osh(2A) = ^ {(A 2 + C 2 )(C 2 + D 2 ) - 2AD + l} . (66) 



6 Physical applications 

The immediate application of this similarity transformation is in one dimen- 
sional crystals requiring repeated applications of the Sp(2) matrix, such as 
multilayer optics [101 E], finite periodic potentials [HI [8] and thus surface 
problems, and laser cavities [33]. For those problems, we are confronted with 
the burden of calculating M N . Now, we have 

M N = (GWG^Y = GW N G-\ (67) 
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where it becomes trivial to calculate W N for all its three cases. 

These days, the diagonalization becomes the central issue when we deal 
with "entanglement" problems in physics, especially when the system cannot 
be diagonalized through a rotation. The system of two entangled harmonic 
oscillators is a case in point [H, [5] . Let us start with the ground-state wave 
function ^ 

ipo(xi,x 2 ) = — / =exp^--{xl + xl)\. (68) 



Here, the variables X\ and x 2 are separable. When they are coupled, the 
wave function takes the form 

ip v (xi, x 2 ) = -j= exp j-i e~ v (x 1 + x 2 ) 2 + e n (x 1 - x 2 f j . (69) 
This wave function can be expanded as (5j [20] , 

ip ri {x 1 ,x 2 ) = ) V (tanh- ) (f>k(xi)(f)k(x 2 ) , (70) 

cosh(r?/2) Y\ VJ 

where (f>k{x) is the normalized harmonic oscillator wave function for the k — th 
excited state. 

This expansion serves as the mathematical basis for squeezed states of 
light in quantum optics [3J, as an illustrative example of Feynman's rest of 
the universe [211 [22], an d more recently as an entangled oscillator state [HE]. 
In order to obtain equation (|69l) from equation (!68|) . we first have to rotate 
the coordinate system with the matrix 



1 /l -1 



(71) 



Then we have to squeeze the new coordinates using the matrix 

(72) 



e-"/ 2 

which is the same as the matrix given in equation f[36|) . Indeed, this system 
was coupled by q squeeze S{rj) followed by a rotation by 45° degrees. 

In addition, this paper examines the question of transition from R{<p) to 
X(x) via N(j) given equation f[35|) . This is known as the stability prob- 
lem in cavity physics [T3] and focal condition in lens optics [12]. While its 
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mathematical aspect was studied in detail in multilayer optics [9J, its physi- 
cal implication is yet to be studied. There could be many other transitional 
problems like this in physics. 

Let us note again that the group Sp(2) is locally isomorphic to the Lorentz 
group applicable to two space dimensions and one time variable. The tran- 
sition from R(4>) to Ni^y) is known as a group contraction, which is similar 
to the transition of the internal space-time symmetries of massive particles 
to that of massless particles in the limit v — > c [17]. It is interesting to 
note that this process can be achieved analytically in terms of the Bargmann 
parameters. This requires further investigation. 

Concluding Remarks 

Two-by-two matrices appear in almost all branches of physics. For matrices, 
the basic question is whether they can be brought to a diagonal form. We 
have addressed this basic question in this paper. We started with the most 
general form of the Sp(2) matrix and its Bargmann decomposition. We 
then showed that it can also be written as a similarity transformation of the 
Wigner matrix, which can be either diagonalized or can be written as a two- 
by-two matrix with unit diagonal elements and one vanishing off-diagonal 
element. 

Also in this paper, we have found that the Sp(2) matrix can be written 
both in the Bargmann form of writing it as a product of three-one parameter 
matrices and also in the Wigner form as a similarity transformation of the 
Wigner matrix. The fact that Sp{2) matrix can be written in these two 
different forms may lead to further interesting results. 
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